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AIAA standards S080 and S081 are applicable for certification of metallic pressure 
vessels (PV) and composite overwrap pressure vessels (COPV), respectively. These 
standards require damage tolerance analyses with a minimum reliable detectible flaw/crack 
and demonstration of safe life four times the service life with these cracks at the worst-case 
location in the PVs and oriented perpendicular to the maximum principal tensile stress. The 
standards require consideration of semi-elliptical surface cracks in the range of aspect ratios 
(crack depth a to half of the surface length c, i.e., (a/c) of 0.2 to 1). NASA-STD-5009 provides 
the minimum reliably detectible standard crack sizes (90/95 probability of detection (POD) 
for several non-destructive evaluation (NDE) methods (eddy current (ET), penetrant (PT), 
radiography (RT) and ultrasonic (UT)) for the two limits of the aspect ratio range required 
by the AIAA standards. This paper tries to answer the questions: can the safe life analysis 
consider only the life for the crack sizes at the two required limits, or endpoints, of the (a/c) 
range for the NDE method used or does the analysis need to consider values within that 
range? What would be an appropriate method to interpolate 90/95 POD crack sizes at 
intermediate (a/c) values? Several procedures to develop combinations of a and c within the 
specified range are explored. A simple linear relationship between a and c is chosen to 
compare the effects of seven different approaches to determine combinations of a; and c; that 
are between the (a/c) endpoints. Two of the seven are selected for evaluation: Approach I, 
the simple linear relationship, and a more conservative option, Approach III. For each of 
these two Approaches, the lives are computed for initial semi-elliptic crack configurations in 
a plate subjected to remote tensile fatigue loading with an R-ratio of 0.1, for an assumed 
material evaluated using NASGRO™ version 8.1. These calculations demonstrate that for 
this loading, using Approach I and the initial detectable crack sizes at the (a/c) endpoints in 
5009 specified for the ET and UT NDE methods, the smallest life is not at the two required 
limits of the (a/c) range, but rather is at an intermediate configuration in the range (a/c) of 
0.4 to 0.6. Similar analyses using both Approach I and III with the initial detectable crack 
size at the (a/c) endpoints in 5009 for PT NDE showed the smallest life may be at an (a/c) 
endpoint or an intermediate (a/c), depending upon which Approach is used. As such, 
analyses that interrogate only the two (a/c) values of 0.2 and 1 may result in unconservative 
life predictions. The standard practice may need to be revised based on these results. 
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* Senior Scientist 
: Consultant, Member ASTM/ASME. 
* NASGRO® 8.1 is a Fracture Mechanics and Fatigue Crack Growth Analysis Software developed and distributed under the 
terms of a Space Act Agreement between NASA Johnson Space Center and Southwest Research Institute. 
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I. Introduction 


IGHT weight metallic pressure vessels (PV) and Composite Overwrap Pressure Vessels (COPVs) are used in 

spacecraft and launch vehicles to carry fuels and other necessary fluids. The fluids carried by these vessels are 
under high pressure and failure of a PV or a COPV will likely lead to loss of mission for robotic missions and loss of 
mission and life for human rated missions. Hence, the structural integrity of the PVs and COPVs is important. 
During the past three decades, considerable progress has been made in designing, operating, and predicting failures 
of these PVs’. 

In this paper, some observations on damage tolerance (DT) of PVs and COPVs are presented. Linear elastic 
fracture mechanics (LEFM) is used to evaluate life and residual strength of the metallic liner. As is customary in 
fracture mechanics, stress analyses are performed for the maximum design pressure (MDP) in the vessel. Regions of 
high stresses and large stress gradients are identified. In the region of the highest stresses, semi-elliptic surface (part 
through) cracks are introduced and oriented perpendicular to the maximum normal tensile principal stress. The 
growth of these cracks is then evaluated for spectrum loading that the PV experiences in its lifetime. While the 
processes that need to be followed are clear from the AIAA standards S 080° for metallic PVs, and S 081° for 
COPVs, implementation details are often fuzzy and are open for interpretation by the users/analysts. This paper 
examines a few areas where implementation issues arise and offers some guidance with regard to LEFM 
calculations, assumed initial detectible crack sizes, fatigue and fracture considerations (i.e., safe-life) and the use of, 
and limitations of, LEFM. 


II. Pressure Vessel Configurations 


Figure 1 shows a cylindrical metallic PV and two COPVs. COPVs have a metallic liner and a composite 
overwrap that is wound over the liner. Figure 2 shows idealized cross section views of a COPV. The pressure 
vessels can be formed by welding parts of metallic segments together as shown in Fig. 3, or other techniques (e.g., 
spin forming). As mentioned previously the PVs need to comply with standards and requirements laid out in AIAA 


S080 and AIAA S081A standards. —_ - 


Noe | 


Hemispherical Domes 


(b) Cylindrical Composite Overwrapped Pressure Vessel (c) Spherical Composite Overwrapped Pressure Vessel 


Figure 1. Example of metallic pressure vessels (PV) and composite overwrapped pressure vessels (COPV). 
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Figure 2. A COPV and its cross section. Figure 3. Welded pressure vessel. 


III. Damage Tolerance Requirements 


The damage tolerance requirements in the standards specify that all areas of the metallic liner except the shear 
region of the boss shall possess a minimum damage tolerance life (safe-life) of 4 lifetimes under the spectrum 
loading the PV experiences during its lifetime. The standard also requires demonstration of residual strength at 
MDP. That is, the maximum stress intensity factor at any point on the crack front of the semi-elliptic crack must be 
less than the condition resulting in fracture of the material, that may be the fracture toughness of the material, Kic, or 
other relevant parameters listed in the later section titled “Other Considerations.” The assumed initial flaw (crack) 
size and shape for both safe-life and residual strength is based on nondestructive evaluation capability (i.e., a crack 
size with an estimated probability of detection of 90% with 95% confidence). The NASA-STD-5009 Table 1* shows 
the NDE capabilities and the corresponding size of the cracks that can be reliably detected with “Standard” NDE 
methods and qualified inspectors. A version of this table is shown in Table 1. The data is derived from historical 
studies that estimated the 90% probability of detection at a 95% confidence level for detecting cracks in flat plates 
by the NDE methods as described in the standard. A crack that is the size estimated by the 90/95 POD for the given 
NDE inspection is presumed to exist in the structure. Different NDE techniques may be used at different locations 
on the vessels. Thus, the size and shape of the initial flaw that is used to demonstrate safe-life can be different at 
different locations. In addition, the safe-life demonstration must consider how the thickness and the material 
properties vary at different locations (e.g., acreage, attachment bosses, welds). Also, localized stresses may include 
significant bending as well as tension stresses in regions of geometric shape changes or curvature, and in addition, 
welding residual stresses may be present at welded joints. All these variables need to be accounted for in the DT 
analyses. For the analyses reported here, the effects of a set of simple remote tension and bending loading examples 
are evaluated and discussed. 

Table 1. NDE Methods and detectability limits. 


Crack Location Part Thickness, t, in. Crack type Crack Dimension, a, in. Crack Dimension, c, in. 
Eddy Current 
Open Surface t<=0.050 Through t 0.050 
t>0.050 PTC 0.020 0.100 
0.050 0.050 
Penetrant NDE 
Open Surface t<=0.050 Through t 0.100 
0.050<t<0.0750 Through t 0.150-t 
T>0.0750 PTC 0.025 0.125 
0.075 0.075 
Radiographic NDE 
Open Surface t<=0.107 PTC 0.7t 0.075 
T>0.107 PTC 0.7t 0.7t 
Embedded 2a=.7t 0.7t 
Ultrasonic 
Open Surface t>=0.100 PTC 0.030 0.150 
0.065 0.065 
Embedded 0.017 0.087 
0.039 0.039 
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IV. Damage-Tolerance Life (Safe-Life) Analysis 


Linear elastic fracture mechanics (LEFM) is generally used to demonstrate damage tolerance life for metallic PV 
and COPVs with elastically responding liners. A part-through crack is assumed to be in the most critical locations 
and is oriented perpendicular to the highest principal stress in the vessel. The flaw (crack) shape (a/c) in the range of 
0.2 to 1.0 are considered (see Fig. 4). 

A state-of-the-art crack growth software package 
such as NASGRO*» is generally used to conduct the 
analysis. As the crack grows under cyclic loading, 
the flaw shape changes are accounted for in the 
analyses. Details used in the LEFM analyses such as 
plane strain fracture toughness (Kj.) and fatigue crack 
growth rates, (da/dN); the loading spectrum and /}———— 2c -=— c-tip 
environments; NDE methods and corresponding 
initial flaw sizes; and a summary of significant 
results are reported in the DT analyses. Both AIAA S080 and S081 standards also allows for the demonstration of 
DT life through testing of coupons or flight representative COPVs. 


a-tip 


Figure 4. Semi-elliptic surface crack. 


Vv. Assumptions 


Damage tolerance safe-life analyses, as is practiced today, has the fundamental assumption that LEFM is valid. 
The underlying LEFM assumptions are: 
1. The material of the vessel away from the crack tip is undergoing elastic deformations. 
2. The plastic region ahead of the crack tip/front is small in relation to the length of the crack. 
3. Deformations outside the small plastic regions such as the remaining ligament between the crack front 
and the far surface of the thickness are presumed to remain elastic. 
4. The material is presumed to have a grain structure such that the integrated effects of the grains in these 
regions acts like a continuum. 

These LEFM assumptions are generally valid when the crack depth is small compared to the thickness. 
However, the limitations of these assumptions (i.e., what is considered small for the ratio of crack depth to 
thickness) have not been established by the standards. Thus, the user is responsible to verify that LEFM assumptions 
are appropriate for the analysis being conducted. The violation of the LEFM assumptions, as will eventually occur 
as the liner thickness is reduced, may lead to unconservative results. This paper presumes conditions of LEFM 
applicability. 


VI. Critical Question 


While the standard appears to be clear on the procedures/processes needed for demonstrating safe-life, various 
questions arise when implementing these requirements. For example, how is the requirement to examine cracks ‘in 
the range of (a/c) = 1.0 and 0.2’ interpreted? Does the stress-intensity factor need to be examined at every point on 
the crack front? Where is the stress-intensity factor maximum along the crack front? Are the conclusions drawn for a 
particular liner region thickness valid for the entire liner? If not, what range are these conclusions valid? The 
remainder of the paper attempts to answer some of these questions. 

Maximum K locations: Where along the crack front of a part through crack (PTC), also called a surface crack, is 
the stress-intensity factor maximum? The maximum stress-intensity factor (K) location varies with aspect ratios 
(such as (a/c) = 0.2, 0.4, 0.6, 0.8, 1.0) and with crack depths (such as (a/t) =0.2, 0.4, 0.6, 0.8)°°. The maximum K 
location also depends upon the applied loading (tension and bending, etc.). From previous work in references 6-8, 
this maximum K location was identified to occur only at either the c-tip or the a-tip for the tension loading 
condition. Figure 5 (a) — (d) presents the normalized values of K for remote tensile loading, S;, in a flat plate for 
various crack depths (a/t); (a/t) = 0.2 represents a shallow crack while (a/t) = 0.8 represents a deep crack. In these 
figures the value of K at the c-tip (@=0) and a-tip (¢=7/2) are plotted. Figure 5 demonstrates that for (a/c) = 0.2 and 
0.4 the maximum value of K occurs at the deepest point, a-tip (¢=7/2) for all crack depths. In contrast, for (a/c) =1 
the maximum value always occurs at the c-tip, (¢=0). For the crack shape corresponding to (a/c) = 0.6, however, the 
behavior is mixed. For (a/f) = 0.2 through 0.6, the maximum value of K occurs at the deepest point (¢=7/2) , while 
for (a/t) = 0.8 the maximum value occurs at the c-tip, (¢=0). This also demonstrates that the maximum values occur 
at either the a-tip or the c-tip and not elsewhere on the crack front. Thus, the DT life must only be interrogated at the 
a-tip and the c-tip. 
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Figure 5. Normalized Stress-intensity factors for various crack shapes at c- and a-tips for remote 
tensile loading, S. 

Loadings: For other loading cases, the behavior depends on the combined loadings. Examples reported here 
include a pure tension, S;, and several combinations of tension S; plus pure bending, Sy. Figure 6 illustrates these 
loadings that cause crack opening at the c-tip, and more complex conditions at the a-tip. In practice, either S; or Sp or 
both may be reversed, causing compression on the crack c-tip. More complex stress fields also may be encountered, 
such as those due to residual stress effects. For these cases, computer analyses must be used that compute the K 
values due to the integrated effect of local stresses acting at the c-tip and a-tip. These K values are then used to 
predict crack growth rates at each crack tip, the incremental crack growth at each tip, and the resulting progressive 
change of crack shape during fatigue crack growth. 


(a) Surface crack in a plate subjected to (b) Remote tension S; and bending S, loading 
remote tension, S, and bending My). on the top and bottom faces of the plate. 
Figure 6. Surface crack in a plate with tension, S,, and bending loadings M). 
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VII. Crack Shapes and Sizes 


NASA-Standard-5009* discusses NDE requirements for fracture critical components. All pressure vessels and 
COPVs are classified as fracture critical’, The 5009 NDE standard lists several NDE methodologies and the 
limitations of their ability to detect part through cracks (surface) and embedded cracks. The NDE methodologies 
addressed here are eddy current (ET), penetrant (PT), radiography (RT), and ultrasonic (UT). Table 1 summarizes 
the 5009 NDE standard for each of these NDE techniques and the associated minimum detectible flaw sizes (based 
on an estimated 90/95 POD) for through the thickness cracks, surface cracks, and embedded cracks. Each of the 
NDE techniques has a different minimum detectible flaw size and each of these techniques is applicable to particular 
thicknesses of metallic PVs and COPV liners. Surface cracks could be on the internal or the external metallic 
surfaces. Embedded cracks at joints could be centrally located through the thickness of the metal, or could be offset 
from the center. The DT analyses need to account for all these cases. 

Table 1 provides the minimum detectible flaw sizes at the extreme aspect ratios, or endpoints in this analysis for 
“Standard NDE” per 5009”. The safe-life requirement is that the most damaging crack configuration be considered. 
None of the standards describe how to interpolate between the NDE detectable crack aspect ratio endpoints, 
(a/c) = 0.2 and 1, to determine the most damaging crack configuration. Further, there is little or no known POD data 
for the various NDE methods at intermediate (a/c) ratio values that can be used to base any interpolation between 
the endpoints. As such, seven approaches to determine initial crack sizes between the (a/c) endpoints are studied. 
These simple approaches are described below as to whether they are conservative relative to each other. Other 
considerations are discussed later. These approaches are illustrated in Fig. 7 for the penetrant NDE method 
detectable crack sizes at the crack aspect ratio endpoints in 5009. 


Approach I:a=mce+r 

This approach assumes that the crack depth (a) is a linear function of the crack length (c) between the two (a/c) 
ratio endpoints, (a/c) = 0.2 and 1. The slope m and the intercept r of the straight line are derived from the crack 
depth and length at the two endpoints. The resulting straight line is more conservative than Approach II but less 
conservative than Approach III (discussed below). This method was originally proposed in the Space Shuttle Orbiter 
Fracture Control Plan'”, 


Approach I: a = m (a/c) +r 

This approach assumes that the crack depth (a) is a linear function of the aspect ratio (a/c) between the two (a/c) 
ratio endpoints, (a/c) = 0.2 and 1. The slope m and the intercept r of the straight line are derived from the crack 
depth and aspect ratio at the two endpoints. The resulting curve has a non-linear relationship between (a) and (c). 
This relationship produces smaller interpolated crack sizes and is therefore nonconservative relative to Approach I. 


Approach II: c = m (a/c) +r 

This approach assumes that the crack length (c) is a linear function of the aspect ratio (a/c) between the two a/c 
endpoints, (a/c) = 0.2 and 1. The slope m and the intercept r of the straight line are derived from the crack depth and 
aspect ratio at the two endpoints. The resulting curve has a non-linear relationship between (a) and (c) and this 
relationship is conservative relative to Approach I because it produces larger interpolated crack sizes than Approach 
I. 


Approach IV: Area=matr 

This approach assumes that the crack area (% za c) is a linear function of the crack depth (a) between the two 
(a/c) endpoints, (a/c) = 0.2 and 1. The slope m and the intercept r of the straight line are derived from calculating the 
area at each endpoint and fitting a line to the area and crack depth. This approach produces smaller interpolated 
crack sizes and is less conservative than Approach I. 


Approach V: Areaa=mct+r 
This approach assumes that the crack area (% zac) is a linear function of the surface crack length (c) between 
the two (a/c) endpoints, (a/c) = 0.2 and 1. The slope m and the intercept r of the straight line are derived from 


° Other, smaller crack sizes can be considered if the NDE techniques and operators have passed a POD evaluation to prove 
capability as “Special NDE” per requirements in section 4.3 in the NASA-STD-5009%. It is recognized that most thin wall flight 
PV are qualified for “Special NDE” to obtain smaller detectability limit crack sizes than are in 5009 Table 1. However, “Special 
NDE” crack sizes are not published, as they are individually established per 5009 requirements, thus cannot be used for the 
example analyses reported in this paper. Therefore, the “Standard NDE” crack sizes are used for the assessments reported here. 
Nevertheless, the analyses reported here are expected to also be applicable for “Special NDE” initial crack sizes, although the 
results may vary somewhat due to the different initial crack sizes. 
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calculating the area at each endpoint and fitting a line to the area and surface length. This approach also produces 
smaller interpolated crack sizes and is less conservative than Approach I. 


Approach VI: Area = constant (based on the area for the NDE endpoint at (a/c) = 1) 

This approach calculates the area at the endpoint at (a/c) = 1 and assumes that the area is constant for all other 
aspect ratios. As discussed later in the paper, this approach may be more conservative than Approach I and does not 
necessarily agree with the value at the (a/c) = 0.2 endpoint. 


Approach VII: Area = constant (based on the area for the NDE endpoint at (a/c) = 0.2) 

This approach first calculates the area at the endpoint at (a/c) = 0.2 and assumes that the area is constant for all 
other aspect ratios. This approach is less conservative than Approach I and does not necessarily agree with the value 
at the (a/c) = 1 endpoint. 


0.08 a/c=1.0 
0.07 
0.06 a-tip 
< as ear 
ae, Mavi 
— 2 
0.03 oe 
=e Approach | _ 
0.02 —*—Approach II wee 
0.01 —«— Approach Ill 
0 
0.00 0.02 0.04 0.06 0.08 0.10 0.12 
c, in. 
(a) Approaches I, II, and III. 
0.08 a/c= 1.0 0.08 vis 
0.07 0.07 
0.06 0.06 
: 0.05 0.05 
© 0.04 < 0.04 
© 0.03 ; 
== Approach | © 0.03 
0.02 —s— Approach IV a/c=0.2 0.02 =e Approach | a/c=0.2 
0.01 —— Approach V 0.01 pass 
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Oe et G06 00S: 02) Ele 0.00 0.04 0.08 0.12 0.16 
C, In. c, in. 
(b) Approaches I, IV, and V. (c) Approaches I , VI, and VII. 


Figure 7. Various crack shape interpolation approaches studied using 5009 detectable crack aspect ratio 
endpoints for penetrant NDE. 


VIII. Discussion of Various Approaches 


The lack of minimum detectible flaw sizes for the standard NDE methods in 5009 for aspect ratios between the 
two specified NDE aspect ratio endpoints (a/c = 0.2 and 1.0) requires an assumption be made when predicting life 
for cracks between these endpoint sizes. The NDE information is given in terms of crack depth, a, and crack half- 
length, c, so the choice was made to assume a linear relationship between these two parameters (Approach I) and 
perform calculations with results as shown in Figs. 8 through 12, 14, and 15. All other approaches use quantities 
such as (a/c) or crack area or another assumption. These approaches may be more or less conservative than 
Approach I as detailed previously and in the following discussion. 
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Figure 8. Eddy current NDE 5009 detectable crack aspect ratio endpoints, interpolated crack 
shapes using Approach I, and lives at each aspect ratio (20 ksi remote tension & plate 
thickness=0.075 in). 
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Figure 9. Eddy current NDE 5009 detectable crack aspect ratio endpoints, interpolated crack 


shapes using Approach I, and lives at each aspect ratio (20 ksi remote tension & plate 
thickness=0.1 in). 


8 
American Institute of Aeronautics and Astronautics 


a-tip 


| ar 


a 


0.065, 0.065 


0.06 


-®— Crack Length: c 


—™-Crack Depth: a 
~~ Life/21,000 


0.150, 0.030 


Crack Depth a, in. 
o 
‘4 
Crack Dimension, in. 


0.00 0.05 0.10 0.15 0.20 


Crack Length c, in. 


0.0 0.2 0.4 0.6 08 1.0 


Crack Aspect Ratio a/c 


Life / 21,000 


(a) Ultrasonic NDE aspect ratio endpoints and (b) Calculated lives for Approach I interpolated 


interpolated shapes using Approach I. initial crack size and aspect ratios. 


Figure 10. Ultrasonic NDE 5009 detectable crack aspect ratio endpoints, interpolated crack 


shapes using Approach I, and lives at each aspect ratio (20 ksi remote tension & plate 
thickness= 0.1 in). 
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Figure 11. Penetrant NDE 5009 detectable crack aspect ratio endpoints, interpolated crack 
shapes, normalized initial stress intensity factors, and resulting lives at aspect ratios (20 ksi 


remote tension & plate thickness=0.1 in). 


9 
American Institute of Aeronautics and Astronautics 


0.5 2.0 


1.8 
0.4 
——Approach Ill: c= m (a/c) +r 1.6 

0.08 0.075, 0.075 —w—Approachl: a=mctr wn 5 
© NDE (a/c) Endpoints ~ 0.3 as) 
cs — Approach Ill: K(a) / S_ £ 
0.07 ice —G=Approach Ill: K(c) / Sy 1.4 < 
: = ——Approach I: K(a) / Sy a 
= 0.06 a —=Approach I: K(c) / Sy | 2 
© =@=Approach Ill Life / 24,000 Ey 
= 0.05 =@=Approach | Life / 27,000 12 > 
a 2 
a £ 
~~ 0.04 a 
6 0.1 wv 
15) 0.125, 0.025 1.0 Y 
0.03 3 

0.07 0.08 0.09 0.10 0.11 0.12 0.13 0.0 0.8 

Creech Gan, 0.0 0.2 0.4 0.6 0.8 1.0 
Crack Aspect Ratio a/c 
(a) Penetrant NDE aspect ratio endpoints and (b) Initial crack tip K(a), K(c) for interpolated 
interpolated shapes using Approaches I and IIL. crack shapes using Approaches I and III, and 


computed lives at each aspect ratio. 


Figure 12. Penetrant NDE 5009 detectable crack aspect ratio endpoints, interpolated crack 
shapes, normalized initial stress intensity factors, and resulting lives at aspect ratios (20 ksi 
remote tension & plate thickness=0.12 in). 
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Figure 14. Comparison of lives for eddy current NDE crack shapes with remote tension (20 
ksi) and tension (15ksi) + bending (5 ksi) loadings. 
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Figure 15. | Comparison of lives for Eddy current NDE crack shapes remote tension (20 ksi) 
and various tension + bending loadings. 


The assumption of linearity between a and c (Approach I) is more conservative than the approaches that assume 
that the area varies linearly between the two endpoints (Approaches IV & V). Thus, Approaches IV and V are not 
recommended. 

The assumption that the surface crack depth is a linear function of the aspect ratio between the two endpoints 
(Approach IT) is less conservative than Approach I. Thus, Approach II is not recommended. Approach HI has the 
advantage that it is more conservative than Approach I, so it is also evaluated with results shown in Figs. 11 and 12 
for the 5009 PT NDE detectable crack aspect ratio endpoints. 

The assumptions that the crack area is constant (Approaches VI and VII) has a logical inconsistency because the 
areas at the two endpoints are generally not the same. Thus, using the constant area of one of the NDE endpoints 
will have a discontinuity at the opposite endpoint. The approach that uses the smaller area (generally the smaller 
aspect ratio: Approach VII) will likely be unconservative relative to Approach I. The approach that uses the larger 
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area (generally the larger aspect ratio: Approach VI) can be either conservative on nonconservative relative to 
Approach I. Further guidance from the NDE community is needed if Approach VI is preferable. 

From this study, it is also noted that it is not clear how to predict the initial crack size for other aspect ratios if 
only one NDE detectable crack aspect ratio endpoint value (e.g., (a/c) = 1.0) is known. The use of an equal area 
assumption would appear to be inconsistent with the detectable crack aspect ratio endpoint sizes with different areas 
found in 5009* 

As noted previously, Approach III is more conservative than Approach I, and hence can be used as the method 
for interpolating between the (a/c) ratio endpoints for all the reported results. However, the results reported for the 
ET and UT NDE methods were performed first using Approach I. Also, Approach I is appealing because it produces 
crack sizes that are on the straight line connecting the two endpoints, i.e., it produces a, c points that are always 
closest to the known a, c, values at the (a/c) ratio endpoints. Both Approach I and III are used to provide a 
comparison of their results for the PT NDE method evaluations described later. 

Discussions with available NDE experts confirmed the 5009 NDE data values represent probability of detection 
data from historical controlled studies as cited in the standard and the data were not derived based on any a priori 
presumptions as to how crack shape (a/c) or area or any other aspects affected the detectable crack sizes. In other 
words, at this time, it appears there is no information suggesting how to interpolate between the (a/c) endpoints. 


IX. Life Calculations 


Fatigue crack growth life predictions were performed assuming that LEFM is applicable. NASGRO*® version 8.1 
was used to perform the life predictions for a surface crack in a 3-inch-wide flat plate of 6061-T6 aluminum. The SC 
30 stress intensity factor solution was selected for this configuration. In the first set of analyses, the plate was 
subjected to a remote tensile stress of 20 ksi with a stress ratio (R = Sinin/Smax Where S is the remote applied stress) of 
0.1. The lives of the crack that are shown in Figs. 8 through 15 are the lives for each of the interpolated initial a and 
c crack size aspect ratios for continuing cycles of the applied loadings up to the point when the crack transitions 
from a PTC (Part Through Crack) to a TC (Through Crack). Also, in these analyses, all easily available NASGRO® 
failure criteria checks (fracture toughness value Kjc, net section stress, plastic limit load, and Option 1 FAD*) were 
used. All analyses predicted failure after the transition to a through crack (TC) for this material due to Option 1 FAD 
calculation; however, the life to the fracture event was not reported because once the PV leaks it has failed to 
perform its function. 

Figures 8 through 12 present initial interpolated crack shapes and the resulting lives for the semi-elliptical cracks 
in a plate subjected to remote tensile stress of S,=20 ksi (see Fig. 6). The four crack NDE methods presented in 
Table | are considered. Approach I that assumes a linear relationship between a and c of the two (a/c) endpoints was 
used to determine intermediate (a/c) crack sizes in all the analyses reported here. However, as noted previously, 
Approach III is more conservative, and results using both Approach III and Approach I appear in Figs. 11 and 12 
that address PT NDE. Note that the Radiographic NDE is an exception that does not need an Approach method for 
interpolating sizes because the 5009 standard specifies that the crack sizes that vary are a function of plate thickness. 
These results are presented in Fig. 13. 

For the other three NDE methods, that are all assessed using Approach I, a typical point j on the straight line 
connecting the a and c values of the endpoints is read off as a; and c;. (Or alternatively, the linear equation for the 
line connecting the a, c endpoints may be used to compute intermediate points.) The corresponding aspect ratio is 
computed as (aj/c;). The combination of aj, c;, and (a//c;) initial crack sizes used in the analysis results are provided in 
plots of a vs c, and their values are also plotted against the left hand side and abscissa of each plot showing crack 
life except Figs. 11 and 12. (Since both Approaches I and III are used for results reported in these two figures, it is 
not practical to show the crack sizes in the plots showing life.) 

In the first procedure for determining crack a and c sizes for analysis as outlined above, round numbers for (a/c), 
for example (a/c) = 0.3 (as opposed to (a/c) = 0.3165), are not easily obtained. If one requires round (a/c) numbers, a 
second procedure may be followed. Simple algebra was used to transform the linear relationship of a and c (Eq. 1) 
into relationships between c and a/c (Eq. 4) and a and a/c (Eq. 5). The a vs c fit can be expressed as 


a=metr (1) 


° FAD=failure assessment diagram 
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the slope m and the intercept 7 can be evaluated using the two end points of (a/c) values of 0.2 and 1, as 


-2 r= Ao2—M* Co2 (2) 
Dividing both sides of the Eq. (1) by c gives 
(a/e)=m +~ 3) 
Rewriting Eq. (3) 
c =r/[(a/c)-m] (4) 


then, once c is known for the selected (a/c), a is easily obtained using 
a =(a/c)*c (5) 


For any given value of (a/c), (for example a/c =0.3), the values of c, can be calculated using Eq. (4) with the 
slope m and the intercept r given by Eq. (2). The corresponding value of a, can be calculated using Eq. (5). The 
values of cz and a, and (a/c), obtained by this procedure plots on the same straight line between the two end points 
as obtained by the first procedure. Example plots of the interpolated a, c values are provided with the estimated 
detectible crack sizes. The crack initial sizes are inputs for the analysis results shown in Figs. 8 through 12. 

The life of the cracked component was calculated using NASGRO® 8.1, as mentioned above, for each selected 
(a/c;) combination. The life curve is plotted using the right hand ordinate in each figure. This crack life is 
normalized by a convenient number that is approximately the minimum life computed for each case. The 
normalizing numbers or “scale factors” cause the minimum life to plot at 1.0 on the plot right size ordinate, but 
actual life is the "scale factor" times the ordinate. 

Figures 8 and 9 for two plate thicknesses show that the minimum life for the initial crack size determined by 
eddy current NDE does NOT occur at the (a/c) range endpoints, but somewhere near the aspect ratio of about 0.6 for 
a 0.075-inch-thick plate and of about 0.5 for a 0.1-inch-thick plate. Figure 10 shows that minimum life for ultrasonic 
NDE initial cracks does NOT occur at the endpoints, but near (a/c) of about 0.4 for the 0.10-inch-thick plate. 

Figure 11(b) illustrates the effects of using Approach II versus Approach I to interpolate crack sizes for aspect 
ratios between the 5009 crack aspect ratio endpoints for penetrant NDE. When Approach I is used, the minimum life 
occurs at the (a/c) = 1.0 endpoint. When Approach III is used, the minimum life occurs at an (a/c) value of 0.8. 
There are two factors in play causing this effect. One factor is the influence of the surface crack (a/c) ratio on stress 
intensity factor (SIF), as shown in Fig. 5. The other factor results from the larger crack a and c dimensions that 
develop from the Approach III interpolated (a/c) ratio endpoints as compared to those using Approach I. The two 
interpolated crack shape curves, shown in Fig. 11(a), demonstrate the larger crack sizes that result from the use of 
Approach HI. Whether these interpolated shapes are realistic or artificial cannot be determined without addition 
NDE data on crack sizes with equal probability of detection at intermediate (a/c) aspect ratios. The integrated effects 
of the (a/c) effects on the SIF, and the differences in crack sizes at each interpolated aspect ratio are reflected by the 
curves showing the SIF, i.e., K(a) and K(c) for both Approach I and III in Fig. 11(b). The SIFs resulting from 
Approach III are larger in the middle of the crack aspect ratio range, and these tend to move the location where the 
life is a minimum towards the center of the crack aspect ratio range. There is also a smaller effect due to plate 
thickness as illustrated by comparing the results shown in Fig. 12 for a 0.12-inch-thick plate with those in Fig. 11 for 
a 0.1-inch-thick plate. 

In summary, the different aspect ratios with the minimum life shown in Figs. 8 to 12 are due to combined effects 
of both the different values for each NDE method of detectable initial crack sizes, a, c, at the (a/c) ratio endpoints as 
specified in 5009 and the method of interpolating intermediate crack sizes between these a/c ratio endpoints. The 
effect of plate thickness differences on the crack aspect ratio with minimum life is small for the analyses reported 
here. However, larger variations of thickness, and much smaller plate thicknesses, may occur in PV, especially in 
flight COPV liners. Plate thicknesses and variations of thickness in those applications may have a much larger effect 
on life and aspect ratio effects. However, those situations are outside the scope of this study. 

Figure 13 presents the crack shape and life calculations for radiographic NDE PTC cracks. Here the crack depth, 
a, is a function of the thickness of the plate. Hence crack initial lengths a and c are plotted against the thickness of 
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the plate in this figure. For ¢ < 0.107 in., surface crack half-length is a constant c = 0.075 in. The crack depth, a, on 
the other hand is a linear function of the thickness (0.7/). For t > 0.107 in., the surface crack half-length c, is also a 
linear function of thickness (0.7t). This combination of crack shapes results in a normalized life that was nearly 
linear up to ¢ = 0.107 in. and decreased sharply for thicker plates. Therefore, the use of estimated minimum 
detectible crack sizes for RT will require a careful evaluation of the influence of plate thickness. 


X. Tension and Bending Loadings 


Computations for combinations of remote tensile and bending loadings for the eddy current (ET) NDE method 
initial detectable cracks are described next. For each of these loading conditions, the same eddy current initial PTC 
sizes using Approach I for interpolation of intermediate crack sizes shown previously in Figs. 8 and 9 are used. 
These PTC sizes appear again in Fig. 14. Four different remote loadings are analyzed including: (A) S; =20 ksi 
tensile loading, (B) S; =15 ksi tensile plus S, = 5 ksi bending loadings, (C) S; =10 ksi tensile plus S, =10 ksi bending 
loadings, and (D) S; = 0 tensile with S, = 20 ksi bending loading. These stress loadings are shown in Fig. 6 (b) as S; 
and Sp». Note that the four cases were selected to have the same maximum stress acting on the outer surface and the 
crack c-tip. These analyses demonstrate the tension loading has a much shorter life than cases consisting partly of 
bending stress. 

Figure 14 presents comparison of normalized lives for the first two cases (A) and (B) for plate thicknesses of 
0.075 in. and 0.10 in. As expected, the pure tensile loading has a lower life for most aspect ratio crack configurations 
for both plate thicknesses considered. The minimum life for the (B) loading occurs at a slightly smaller aspect ratio 
for the two plate thicknesses as compared to the (A) loading minimum locations -- that were previously shown in 
Figs. 8 and 9 and noted to be about 0.6 for the 0.075-inch-thick plate and about 0.5 for the 0.1-inch-thick plate. The 
change from all tension to tension plus a small bending loading shown in Fig. 14 has not produced a significant shift 
in the crack aspect ratio where the minimum life occurs. (Note the “scale factor” causes the Tension only 20 ksi 
loading (labeled “Ten 20”) to appear to have a larger ordinate (life) than the loading with 15 ksi tension and 5 ksi 
bending (labeled “Ten 15 & Bend 5.”). In actuality, the 20 ksi tension result has a much smaller life when the plot 
life ordinate is correctly interpreted using the “scale factors” shown in the plot legend — these factors are the 
minimum life for each analysis.) 

Figure 15 presents lives as a function of aspect ratio of the cracks for the same eddy current initial PTC sizes, for 
the four loading cases listed above. As expected, the lives of all the cracks becomes larger as the tension component 
decreases and the bending component increases, even though the total stress at the crack c-tip is always the same. 
The smallest value of life appears to occur at an (a/c) value of about 0.6 for the 0.075-inch-thick plate and 0.4 to 0.5 
for the 0.10-inch-thick plate. The crack aspect ratio where the minimum life occurs has not changed significantly 
with the increasing bending loading cases studied here. However, for loading cases where the bending is fully 
reversed, or where complex stress gradients occur, it is expected the (a/c) value where minimum life occurs could be 
strongly affected. 


XI. Other Considerations 


The standards S 080-1998 and S 081-2006 have different requirements when safe-life (of 4 service life times) 
and leak before burst (LBB) needs to be demonstrated. In these calculations, the applicable fracture toughness (e.g., 
Kic, Kic, Kie, Kc) of the material should be considered. The definitions of various toughness values are given in 
NASA-STD-5019A” and presented in the box below. The values of the toughness used in the analyses need to be 
justified and shown to be valid for the materials and thicknesses used in the vessels and the crack configuration 
being evaluated. 

The manufacturing and fabrication of the vessels generally leaves residual stresses that may influence the 
behavior of cracks. These residual stresses have complex through-the-thickness distributions and may vary from 
location to location on the PV. Often, the presences of residual stresses are suspected, but the magnitudes and 
distributions are unknown. Methods that bound the influence of residual stresses on DT may determine if 
quantitative characterization of residual stresses is necessary. 
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Various Toughnesses Used in the DT Analyses 


K,: Plane stress fracture toughness. 

The value of stress intensity factor K at the tangency between a crack extension resistance curve 
(R-curve) and the configuration-dependent applied K curve (ASTM E1823, Standard Terminology 
Relating to Fatigue and Fracture Testing). This crack extension occurs under conditions that do not 
approach crack-tip plane strain. The R-curve and K, vary with the material, specimen size, and 
thickness. Ke is used in NASGRO® to represent fracture toughness as a function of thickness for 
use in crack growth calculations. 


K,.: Plane strain fracture toughness. 

The crack extension resistance under conditions of crack-tip plane strain in Mode I for slow rates 
of loading under predominantly linear-elastic conditions and negligible plastic-zone adjustment 
that is measured by satisfying a standardized procedure with validity requirements (ASTM E399, 
Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness K;, of Metallic 
Materials). Another quantity, Ky, defined for conditions with limited plasticity from Jj, may also 
be useful (ASTM E1820). 


Kj: Effective fracture toughness for a surface or elliptically shaped crack. 

The toughness is based on residual strength and the original crack dimensions. This parameter is 
meaningful only when crack-tip plastic zones are small and stable crack growth before failure is 
generally absent (ASTM E740/E740M, Standard Practice for Fracture Testing with Surface- Crack 
Tension Specimens, main body and section X1.2). For conditions with plastic effects and well- 
defined crack-tip stress fields with fracture controlled by crack initiation, an approach involving 
constraint may be applicable (ASTM E2899). Testing of flaws in specimens representative of the 
structure is needed to determine damage tolerance for plasticity conditions when crack-tip stress 
fields collapse. Ky. is used in NASGRO” for analyses of crack growth. 


Keac: 

The largest value of the stress intensity factor at which crack growth is not observed for a pre- 
cracked through-crack specimen of specified material, environment, and thickness that is tested for 
a significant duration in accordance with ASTM E1681. 


Kigac: 

The largest value of the stress intensity factor at which crack growth is not observed for a pre- 
cracked through-crack specimen of specified material, environment, and thickness that is sufficient 
to meet requirements for plane strain and is tested for a significant duration in accordance with 
ASTM E1681. 


Kisce? 
Keac is often denoted as K {gcc in the literature. 


XII. Recommendations 


The results presented in this paper illustrate how the interpretation of the implementation details of the AIAA 
S080 and S081 standards can influence the DT life predictions of COPVs. The requirement that safe-life be 
demonstrated for the minimum detectible flaw sizes for specific NDE methods (as assessed by estimated 90/95 
POD) requires an interpretation of the initial part through crack size between the aspect ratio values (endpoints) 
listed in the NASA STD 5009 (generally at an aspect ratio of a/c = 0.2 and 1). The choice of interpolation method 
for the part through crack sizes between the listed NDE endpoints can influence the analytically predicted life. The 
following recommendations are made based on the findings of this paper: 

e Table 1 of NASA STD 5009 that lists the minimum detectible flaw sizes of various NDE techniques 
should be updated to include a minimum of three aspect ratios for each technique. It is recognized that 
this is a significant undertaking, but it is necessary to eliminate the ambiguity present in the current 
situation. 

e In the interim, i.e. before Table 1 of NASA STD 5009 can be updated, the community should adopt a 
standard approach for the interpolation of part through crack sizes from the listed NDE endpoints. There 
are two possible choices: Approach I or III as discussed in the paper. Without any other information, 


15 
American Institute of Aeronautics and Astronautics 


Approach III, as the most conservative, is a logical choice to be the standard method for interpolation 
between the listed crack length (c) and depth (a) values. 


XIII. Concluding Remarks 


Some observations on damage tolerance (DT) analyses of pressure vessels (PVs) and composite overwrap 
pressure vessels (COPVs) are discussed. The two standards that are applicable are AIAA S 080 for metallic PVs and 
S 081 for COPVs. 

As is customary in any DT analyses, a flaw that is the estimated minimum reliably detectible crack is assumed to 
exist in the pressure vessel and the crack is oriented perpendicular to the maximum tensile principal stress. The DT 
analysis then proceeds to demonstrate safe-life (four lifetimes of the component) to satisfy the standard 
requirements. These analyses should use available failure criteria to ensure the results are the minimum life. The 
minimum detectible crack sizes are given in the NASA-STD-5009 for various NDE methods. In this paper, eddy 
current, penetrant, radiographic, and ultrasonic methods were considered. NASA-STD-5009 specifies a semi-elliptic 
surface crack with aspect ratios (a/c) — a is the crack depth and c is the semi-surface length — of 0.2 and 1 for the 
eddy current, penetrant, and ultrasonic NDE methods. No intermediate (a/c) NDE crack size data are provided in the 
standard. Crack size is specified in this standard as a function of thickness for the radiographic method 

While the AIAA standards appear to be clear on the procedures/processes needed for demonstrating safe-life, 
various questions arise from these requirements. For example, how is the requirement ‘in the range of (a/c) = 1.0 
and 0.2’ interpreted? How many crack aspect ratios (a/c) need to be considered? Does the stress-intensity factor 
need to be examined at every point on the crack front? Where is the stress-intensity factor maximum along the crack 
front? Are the conclusions drawn for a particular liner thickness valid for others? In this paper, some of these 
requirements are highlighted and some clarifications are offered. 

The maximum stress-intensity factor location occurs only at the a-tip or the c-tip and not elsewhere along the 
crack front. Several approaches were examined for initial crack sizes between the a/c ratio endpoints as specified for 
the eddy current, penetrant, and ultrasonic NDE methods. The approaches included: linear relationships between a 
and c (Approach J); a and (a/c); c and (a/c) (Approach III); a and the area of the crack; c and the area of the crack; 
constant area determined by (a/c) = 1; and constant area determined by (a/c) = 0.2. Approach I that assumes a linear 
relationship between a and c at the two endpoints is a simple, historical method, but Approach III is more 
conservative. If no other detectable crack size information is available to guide the interpolation, Approach III is the 
logical choice at this time. However, the ambiguity that results from only the two existing (a/c) endpoints could be 
eliminated if a third (a/c) point is established for each NDE method. Also, in that event, additional interpolation 
methods may be considered. 

Using the cited approaches, interpolated crack shapes were constructed between the 5009 aspect ratio endpoints 
for the eddy current, penetrant, and ultrasonic NDE methods and the lives of each of these crack configurations in a 
plate subjected to remote tensile loading was evaluated using NASGRO® 8.1 software. These calculations using 
Approach I demonstrate that the eddy current and ultrasonic NDE methods 5009 detectable crack sizes cause the 
smallest life to not occur at the two (a/c) endpoints, but at an intermediate configuration in the range (a/c) of 0.4 to 
0.6. The penetrant NDE method has smallest life at (a/c) of 1.0 when crack size aspect ratio is interpolated using 
Approach I. However, the result shifts to an intermediate (a/c) of 0.7 to 0.8 when crack size aspect ratio is 
interpolated using Approach III. The effect of plate thickness differences on the crack aspect ratio with minimum 
life is small for the analyses reported here. However, larger variations of thickness, and much smaller plate 
thicknesses, may occur in PV, especially in flight COPV liners. In general, the crack size aspect ratio where 
minimum life occurs is expected to depend on the detectable crack size aspect ratio endpoints, the interpolation 
method, and the loading conditions. As such, analyses that interrogate only the two endpoints of (a/c) =0.2 and 1 
likely result in unconservative life predictions. The standard practice needs to be reviewed based on the results in 
this paper. 
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